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THE SMALL-DISTURBANCE METHOD FOR FLOW OF A COMPRESSIBLE 
FLUID WITH VELOCITY POTENTIAL AND STREAM FUNCTION 
AS INDEPENDENT VARIABLES 


By Carl Kaplan 
SUMMARY 


The equations of two-dimensional compressible flow are treated 
according to the Prandtl-Busemann small-disturbance method. In contrast 
to the usual procedure, the independent variables are the compressible 
velocity potential and stream function and the dependent variables are 
the rectangular Cartesian coordinates in the plane of flow. The six 
first-order differential equations corresponding to the first three iter- 
ation steps are put into complex-vector form. The particular integrals 
of the resulting set of three equations are then directly obtained. As 
an example, the general results of the analysis are applied to the case 
of subsonic compressible flow past a sinusoidal wall of small amplitude. 


INTRODUCTION 


The problem of the integration of the equations of compressible flow 
past a prescribed solid boundary has been treated in numerous papers and 
in several diverse ways. With the exception of the hodograph- 
transformation method, the equations are usually solved by methods of 
successive approximations. One of these, initiated by Janzen and 
Rayleigh, starts from the incompressible-flow solution and develops the 
compressibility effects in a series of powers of the undisturbed stream 
Mach number. This method is restricted to the subsonic range, since the 
differential equations of the process are always of the elliptic type. 

In general, this method yields the best results for flows past thick 
shapes for which the critical stream Mach number is much less than unity. 
In contrast, the Prandtl-Busemann small-disturbance method of iteration 
in terms of a small parameter starts from the undisturbed flow, the first 
step being the well-known Prandtl linearized solution. This method is 
best suited to thin profiles for which the critical stream Mach number 

is nearly equal to unity and, from the beginning, yields a good approxi- 
mation to the desired rigorous solution. Both the Janzen-Rayleigh and 
the Prandtl-Busemann methods are applicable not only to plane but also 
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to axisymmetric subsonic flows. The Prandtl-Busemann procedure has also 
been utilized for the calculation of plane and axisymmetric supersonic 
flows. A good discussion of these two iterative procedures, including 
a number of basic references, is given in reference 1. 


The author has recently considered the problem of plane incompress- 
ible flow past an infinitely long, sinusoidal wall of arbitrary ampli- 
tude (ref. 2). It was found that this problem could not be treated in 
the physical-flow plane but had to be transferred to the plane of veloc- 
ity potential and stream function. In this plane, the problem was not 
only successfully treated by tbe small-disturbance method but also, the 
solution was rigorously expressed in the form of & nonlinear integral 
equation. This unexpected result has prompted the present investigation 
of the form of the Prandtl-Busemann small-disturbance equations with the 
velocity potential and the stream function as independent variables. 

The equation of continuity and the condition of irrotationality form the 
basis of the procedure. Each iteration step consists then of a pair of 
first-order differential equations of the Cauchy-Riemann type rather than 
a single second-order differential equation of the Laplace or Poisson 
type. These pairs of equations, moreover, can be put into complex- 
vector forms and their particular integrals obtained in a straightforward 
manner. The present paper contains the results of this investigation 
inclusive of the third-order step. The example of subsonic flow past a 
Sinusoidal wall of small amplitude is utilized to illustrate the general 
results of the analysis. No counterpart of the integral-equation approach 
of reference 1 could be found. 


ANALYSIS 


The requirements of continuity and irrotationality lead to the 


relations 
D 
u = d... = a y 
(1) 
Po 
y = = =- 2 y 
r i 
where 
u,v velocity components in direction of x- and y-axes, 


respectively 
X,y rectangular Cartesian coordinates in flow plane 
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D density of fluid 

Da density of fluid in undisturbed flow 
g veLocity potential 

Y stream function 


The subscripts denote differentiation with respect to the designated 
variables. The quantities x, y, Ø, wv, wu, and v are nondimen- 
sional with a characteristic length associated with the solid boundary 
as unit of length and the speed U of the undisturbed stream as unit 
of velocity. 


From Bernoulli's theorem there follows the relation 


l 
E - | - 2-5 Mola? - :) m (2) 
where 
M. Mach number of undisturbed stream 
q nondimensional speed of fluid 
y ratio of specific heats at constant pressure and 


constant volume 


The introduction of Ø and Y as independent variables leads to 
the following relations in place of equations (1): 


Poo 
Ty 77g ^j 

(3) 
4 =- -2y 


For the purpose of deriving the iteration equations of the small- 
disturbance method, let 


En 2 
x =P + ex] + exo + exo... 


(+) 
y = Y + ey, + Yo + yz +. e 


——— — — M —— nnn ELE 
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Where xX, and yp (n-21,2,5...) are functions of the independ- 


ent variables Gg and Y and where c is a suitable small parameter 
Which serves to regulate the course of the iteration process. 


Now, 
2 2. 
o 07 +N 


q^ = 5 (5) 
gry - xg) 


and 
Xj = 1+ exyg + Log + Cag + Ew ae 


Xy = ey + o + sy Fs... 
(6) 
Jd = EY1g + e yag + Tag d 3x 


Yy = 1 + yay + EY oy + € ay + de u 
Substituting equations (6) into equation (5) yields 


q^ - 1 = -eexid + efg + eXyyX1d T Xy. - 2g + 2c? (2x4? = 
RE > ANI t ?ng*og + *1yYog + Saig + 


"Cay 7 gy 7 Xu Yay g xg) anal Tm 
Also, from equation (2), 


2 =% 5 — 2 6 


5 
A8 » -1) RES (8) 
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By means of equations (5), (6), (7), and (8), equations (3) yield the 
following sets of relations when the coefficients of the various powers 
of e are equated to zero: l 


pg = Fry | 


(9) 
Xy = Vig 
B^ x. = Yoy + lw Yid. _ |^ + (y + mf! 
(10) 
2 
Xoy = Yag + Ma “Pig 
p xag = Vay ES al 19 og - d + Les + LIC v.) xig e 
2 2 1,4 2 
|e + (y + 1), xs gto + E + (y + 13M," - 
2Q EP + a $ se (11) 


211 3 
zay = -yag + Moo p vag + Xoplig * P2 - 


a 4 "a gp. = E xag 
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By defining 
Xp + $ Yn = 7n Xn - + ya = Za 
Ø + ipp =w d- ipy = 


equations (9), (10), and (11) can be expressed in the following complex- 
vector forms: 


N (12) 


n 
TT + 1 M. 2 — L 2 + 1 Mo TN 
tom = HEURE s + aie) bez) 17 


eae es a jus 1,2 +1 Ale 2 — 2 
257 = - bM. Tay (Zow + Zan) - g M (1 + g + L ma ass - 215 J- 


n 
= M ^p? (215 = Zam) 2 so E 1a- (enr + zs = + Zog + zou + 


2 Y 
= lx e$ M, 5 Mo, 
xri dido ae ere 


h 
au + s (zi " zs) (14) 
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Corresponding to these equations are the following second-order partial- 
differential equations: 


Zy = 0 (15) 

i 

1 M, 
Zoya = - PE 2 nud (16) 

y 
L Ma 2 
ewe =~ NE + Hag - 

EN p^ + (y + Ds 214g + X1j22gd (17) 


Note that, when Ø and y are the independent variables, the basic 
iteration equations (12), (13), and (14) for x and y are of the 

first order. In contrast, when x and y are the independent variables, 
the basic iteration equations for Ø or y are of the second order and 
are similar to equations (15), (16), and (17). 


CALCULATION OF PARTICULAR INTEGRALS 


Equations (12), (13), and (14) can be integrated without difficulty. 
The general solutions are then obtained by the addition of arbitrary func- 
tions of w to the particular integrals. These arbitrary functions are 
of the nature of constants of integration and are chosen in such a way 


that the boundary conditions of the problem under consideration are 
satisfied. 


Both the particular and the general solution of equation (12) there- 


fore is simply that 2, is a function of w only. A particular integral 
of equation (13) can then be written as follows: 


= 2 > 
25 = = ^8 gi — W)Ziy + (zi + je. == 


2 
iw i + Lil E 3 | Zy dv (18) 


e —— — A 
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By substituting this expression for Z» into equation (1+), the fol- 
lowing particular integral can be obtained: 


Y 2 
M M 
2 = Em. (1+ e?) [aan 133 a ps zit m + + 


2 
S (v 2 naase B^ 4 (y + 1) Il 
p 


h} 
TIE zn ein Be 


"M 


2 m= 
(y +1) ——[z zu. os De Bub ydq ee 
g? g? 


TI 


4 2 
M. 2 M 
y (7+1) Sl fs — dW + >. B? + E (7+1) S fa + 


2 
y + 15M : a 


Eo Ys asas (19) 
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The analysis thus far has not been restricted and applies equally 
well to both subsonic and supersonic flows. Actually, of course, Xy 


and yj as functions of f and y are radically different for the 


two types of flow and Be = Me - 1 in the case of supersonic flow. In 


order to illustrate the application of the general results to a partic- 
ular flow problem, the following section contains a treatment of the 
problem of subsonic compressible flow past a sinusoidal wall of small 
amplitude. 


SUBSONIC FLOW PAST A SINUSOIDAL WALL 


The equation of the solid boundary is 


y = € COS x (20) 
and the flow is assumed to be in the direction of the positive x-axis. 
The undisturbed stream velocity U and the wave length ^ in radians 
are utilized as units of velocity and length, respectively. The param- 
eter of the problem is e, the amplitude of the wavy wall in radians. 


The appropriate solution of equation (12), vanishing for Y = œ, is 
simply 


Where A is & real quantity to be determined by the boundary condition 
at the surface, namely, UU =O and y= e cos x. Then 


x= - cA sin g 
and 


cos X = BA cos g 
To the zeroeth order of e, cos x cos Ø and, therefore, A = i or 


x=$- Se Moin g 


(21) 


Dm 


y = 5 V + ce cos g 
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With 2, = 7 elW, the appropriate general expression for Zo can 
be obtained fram equation (18) and is 


i _ 
Zo = 1 7 5 E [or + al - | s 


where B and C are real quantities to be determined by the boundary 
condition at the surface, that is, y =0 with y = e cos x. Therefore, 


2 2 
x=$- ¿sing - (5) guess Lii + p^B|sin 2¢ 


p B 
(22a) 
and 
| l2 Me 
cos x = cos Ø + £ - = Mo cos 2d + (y + 1) (1 - 2 cos 2d + 
8| 8 g^ 4 
B^B cos 26 + BC (225) 


From equation (22a), to the first order of .fp, 


cos x = cos Ø + 5 S (1 - cos ed) 
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Hence, comparison with the right-hand side of equation (22b) shows that 


d 52 g2 
and 
2 1 1 MA 
2 8 


ho ge 
(23a) 
and 
l di M ? 
=y==yt+FeP¥oos g + Ł/E L+ ttie .&-2BW|y +1 Me , 
P B B 4 g2 Y g^ 


E cos 2g (23b) 


- __—_  __Q “Pr —_ — ———— ——— -— — 
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Again, with 2) = - gu equation (19) leads to the general expres- 


sion for Z3: 


2 n 2 - 
| i M. 21.-31W y + 1 Mo 2 y +1 Mo || i(w-2w) 
25 = B (1 + pS)e +i = os M, E A E | 2 


£ 2 y 
1+ pe "ly Hor pla (y 1) 2 + 2 (y + 1) Fa | pedi - 
2 3 ge 8 p? 


2 4 = 
patra tó E eaten : e. p^ + 
go 32 | 


2 À B 
oen S pen fe aem - I 


2 ii E 
2 (y.1)- e ITI "el + Opp) eY pe _ 
9 E ko 


(1 + opy)el(w-2")| tap + inet (21) 


where D and E are real quantities to be determined by the boundary 
condition at the surface, that is, y=0 with y = e cos x. Therefore, 


2 y 
x= $ -$ sing + Tà) E Porph en genns 
P 
(25a) 
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and 


e 2 
cos X = COS f + 5$ (L - cos 2g) + ($) - 5 Me (1 + p“)cos 3d + 


h. 2 4 
z+ UM J| y zg. Mo 1 Mo 
2 2 E 2p 3(y + 1) 2 ae (y + 1) d Ø + 


2 
age 840 e 36 $ + B?D cos Ø + B%E cos 34 
p 
(25b) 


Now, from equation (25a), to the second order of e/B, 


n 
2 
cos x E 1 - &(5) ED cos $ - = cos 2g + 


Comparison of this expression with the right-hand side of equation (25b) 
shows that 


4 2 h 
l gl 2325 2 i Mo _ 1 Meo 
re? * ^5» 32 E + 5(y + 1) 32 5 (y +1) y 


a 
rr A MÀ a A A 


meee o ——— —— ——X—9 ee 
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Tabet phy rt ter _2p2,2 ME 
a *a&^ - a Pt "Xd a tar 


Thus, to the third order of eff, 


2 à 
x= - $ e-BYsin d + (s) : + pe - res = (egy + i e^?P' sin 2g + 


5 
€9[11,1,2.2*1M. 
(5) (Jk * 8 "3 pes * 3 42) 3- 
i 
1 Mo | 4 Qn BY (242 2. 1 gh|.-5B 
bon Be e P'sin ¢ "T TENE sin 3¢ + 
T. 42 
LE ur Me - g? 4 (y & 1) o EG en te n sin jj - 
p? 


2 
bie don ben lamg- 


2 h 
dern dont Melon sa sg - 
g^ p? 


l 
(y + 1) = (py) sin 3g (26a) 
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and 
2 "Ru l 
Ly=Lly+£e Ptos g+ Lfe) Ji z1 ¿ebro Mt 
B B B 2\B g? h p? 


4 2 4 
cale w-3B- OY EL y ae cos Ø + 
= 3i ga 2 p^ 


us 
ÓN -|h + 382 + 


Lyi pe Peo 36 + 27 = : 
h 8 3 8 


2 hc 4 
M 1 MES M 
fan it EA J Bosne woos $ 


2 L 
ven prb ned el cos 3G + 


h 
(y + 1) e (8y )?cos zø (26b) 
B 


— ——— — — A — en re —— EA — -o——— o —À Á— o aaa a^, I Em NE e "EI o lr ur s aru gm ms ur 


-ær —— A A 
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For constant values of the stream function y Z O, equations (26) com- 
prise the parametric equations of the streamlines. The corresponding 
equations for incompressible flow are obtained by putting M_=0 or 


B = 1; that is, 


x= @ - ce" sin Ø + e esin 2d + f(e tsin Gg - eN sin ag ) +.. 
(27a) 
and 
y = Y + ce" cos G+ > e - e^?" cos ag) - 
Z e? (ecos G - e"? cos ag ) : PE" (27b) 


As an exemple of the comparison of streamlines in incompressible 
and subsonic compressible flows, consider the case of e = 0.64, 
M=0.60 or f=0.80, and y = 2.5. The choice of 0.64 for e would 
normally be too large if flow properties were to be calculated at the 
boundary itself with the aid of equations (26). Because of the presence 
of the exponential terms in these equations, however, a rather distant 
streamline such as = 2.5 can be accurately determined. Thus, equa- 
tions (26) become 


x =  - 0.0495 sin Ø + 0.0030 sin 29 - 0.0008 sin 34+... 
and 
y = 2.7865 + 0.0772 cos Ø - 0.0036 cos 24 + 0.0005 cos 59 +... 


Equations (27), for a neighboring incompressible streamline given by 
y = 2.6, become 


x = Ø - 0.0402 sin $+ 0.0011 sin 26+. . 
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and 
= 2.8048 + 0.0402 cos Ø - 0.0011 cos 29 +... 


A comparison of these two sets of equations shows at a glance the well- 
known general result that disturbances caused by the presence of a solid 


boundary die out more slowly for compressible flow than for incompress- 
ible flow. 


CONCLUDING REMARKS 


The use of the velocity potential g and the stream function y 
as independent variables is not new in aeronautical literature, A. Thom 
and his students (particularly L. C. Woods) at the Oxford Engineering 
Laboratory have for years utilized incompressible ø and Y as inde- 


pendent variables with the velocity vectór components log A and 0 


(hodograph plane) as dependent variables. Numerical relaxation methods 
were especially developed for the ¢,-plane and employed in the solu- 
tion of a number of diverse problems in both incompressible and campres- 
sible flows, including mixed subsonic and supersonic flow. In the 
present paper, the investigation is completely analytical and is confined 
to the physical-flow plane. In fact, the compressible velocity potential 
and stream function are thought of as presenting a set of orthogonal 
curvilinear coordinates associated with the shape of the solid boundary 
in the flow plane. The solution of a given flow problem is then expressed 
in the form of a pair of equations transforming the dependent rectangular 
Cartesian coordinates x and y into the independent orthogonal curvi- 
linear coordinates Ø and y. 


Langley Aeronautical Laboratory, 
National Advisory Committee for ies, 
Langley Field, Va., May 17, 1954. 
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